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NOTE ON EINSTEIN'S EQUATION OF AN ORBIT. 

By F. Morley. 

The differential equation of an orbit is given in Eddington's "Report on 
the Relativity Theory of Gravitation"* as 

(1) T A IT = 2m / r 3 - 1/r 2 + 2m/ Wr - (1 - c 2 )/h 2 , 

r and <p being the polar coordinates of the planet. Forsyth f has discussed 
the equation in connection with elliptic functions. The matter is of so 
great interest that a second elliptic handling seems worth publication. 
We have, from Eddington, 

c 2 = 1 — m/a, h 2 = ma(l — e 2 ), 

where a is the "major semi-axis" and e the eccentricity; therefore 

1 _ c 2 = m 2 (l - e 2 )/h 2 . 

The number m?/h 2 is, in Eddington's units, about 10~ 8 . Call it a, so that 

1 - c 2 = a(l - e 2 ). 

In the calculations o? will be neglected when lower powers are present. 
Writing x = m/r, the equation becomes 

(2) (^ Y = 23? - x 2 + 2ax - a(l - e 2 ). 

The roots of the cubic on the right are \ , 0, when a = 0. Hence writing 
x = f + hoc, k 2 a, hza and determining the coefficients h, the roots are 

Xi = \ + 2a, X2 = a(l — e), a; 3 = a(l + e). 

To bring equation (2) into Weierstrass's form 

we write d<p = -\l2dv, x = p + 1/6. The invariants gi and g% are for a 
quartic with coefficients a, b, c, d, e 

ae - 4bd + 3c 2 



* 2d Edition (1920), p. 50. Einstein, Berlin Sitzungsberichte, 1915, p. 831. 
t Proc. Royal Society, series A, vol. 97, April, 1920. 
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Here a, b, c, d, e are respectively 0, 1, — •§, a, — 2o?(l — e 2 ), so that, 
accurately, 

(3) <72 = i(l - 12a), 

(4) gz = sVD - 18a + 54a 2 (l - e 2 )]. 
Hence the discriminant A, or gl — 27 gl, is 

(5) A = 4e 2 a 2 - 64a 3 . 

The roots being real, the network or lattice which gives rise to the 
function p is rectangular. The essence of the problem is the determination 
of a period-rectangle. Let 2«i be the real positive period, 2Xiwi the imag- 
inary period, and write as usual 

a = e~" x . 

The number X is positive, and it may be called the shape of the rectangle. 
The invariants g%, g$, are connected with the periods 2a>i and 2w 2 by 
the equations 

g2 = 60 £' (2m lWl + 2m 2 a> 2 ) 4 ' 9 * = 14 ° ^' (2m lWl + 2m 2 co 2 ) 6 ' 

summed for all integer pairs mi, m 2 , the pair 0, excluded. For calculation 
the double sums are expressed as single sums or products; and the appro- 
priate formula? are given in works on the elliptic functions.* But as the 
proofs of the full formula? are necessarily complicated, I shall interpolate a 
proof of the approximate formula? of a kind that is at once intelligible. 
We have, n being any integer, 

A 1 ir 2 



Differentiating, 



Differentiating again, 



(x + n) 2 sin 2 irx ' 

1 7T 3 COS 1TX 



~£ (x + n) s sin 3 irx 



1 cos 2 irx + i sin 2 irx 



— -7T-4 



|_ sin 4 - 



^(x + ri)* W sirfirx " L sm4ira; 3 sin 2 ira;. 

This enables us to sum the series for g it when m 2 is fixed. When m 2 = 0, 
* They may be found in Harkness and Morley, "Theory of Functions," p. 321, formulas 
75 and 76, and p. 324, formula 82. 
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we get the first tier 

(2«i) 4 L ^2 4 ^3 4 ^ J 2cof 90 12VcoJ" 
When m 2 = 1, we get the second tier 

V. 1 = _J_V___i 

X" (2co 2 + 2m lW i) 4 (2coi) 4 ^ (tX + mi) 4 

= (jlVM *\ 

\ 2wi / \ sin 4 -Kih 3 sin 2 7rzX / 
Now 

sin ^ = ( e " e ) = ( ? - |)/2£. 
Hence the tier m 2 = 1 gives 

6 ofjLVr W + _J1 2 1 

or, neglecting <f, 

W (=)V 

For the tier m 2 = — 1 we get the same result. For the tiers m 2 s ± 2 
we have to replace q by q 2 and so on. Hence the formula 

(6) (^)^2 = 1^ + 20^+ ••-. 
A similar argument, with further differentiation, gives 

(7) (f) 6 <73=^-f<? 2 +--- 
and from these 

(8) (^) 12 A = <f+.... 

From the formulae (6) and (8) we have approximately 

A/fS = 12 3 g 2 

and on the other hand from (3) and (5) we have 

A/gl = 4-27-eV. 
Hence approximately 

q 2 = e 2 a 2 /A 2 , 

q = ea/4. 
For Mercury the eccentricity e is approximately 1/5, so that 

q = 10-720. 
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Calculating irX = log 1/q = 9 log 10 + log 2, the shape X of the rectangle 
of the periods turns out to be nearly 7. As v travels round the rectangle 
of half periods, p takes all real values, so that the radius vector 

r = m/x = m/(p + ^) 

also takes all real values. For the actual orbit we take v on the upper 
side of the rectangle, from v = wi + w 2 , p = a(l + e), corresponding to 
perihelion, to v = co 2 , p = ail — e) or aphelion. 

Finally, the corresponding change in <p is <p = "V2wi, where 



(?)'-* 



for q 2 , being of order a 2 , is here negligible. But from (3) 

0, = i (1 - 12a). 
Hence 

i(^) 4 .i(l-12a) = ^. 

1/(1 - 12a) = 1 + 12a, 

- = 1 + 3a. 

So that in passing from perihelion to aphelion the angle passed through is 

not ir, but 

tt(1 + 3m 2 /A 2 ). 



(5)'- 



